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We study pion-nucleon scattering with a chiral Lagrangian of pions, nucleons, and A-isobars to 
order Q 3 , where Q is a generic small momentum. We compare the results from heavy baryon chiral 
perturbation theory with those from the infrared regularization scheme. While the former provides 
a reasonable fit to the data out to a pion cm. kinetic energy of 100 MeV, the latter is only able 
to fit up to 40 MeV and even then the parameters obtained are unreasonable. Difficulties with the 
infrared scheme in the u-channel are discussed. 



^ ", PACS number(s): 11.30.Rd, 12.39.Fe, 13.75.Gx,13.85.Dz 

O 
O 

(N 
o 

Q 



X 



I. INTRODUCTION 



Pion-nucleon scattering is a fundamental process that one would like to describe using the low-energy realization 
of quantum chromodynamics, namely chiral perturbation theory [1,2]. This is an attractive approach because it not 
only embodies chiral symmetry, which is fundamental to low-energy physics, but also offers a systematic expansion 
in powers of the momentum. Further it ensures unitarity order by order. Gasser and Leutwyler [2] have shown that 
. chiral perturbation theory works nicely for mesons and high accuracy can now be achieved for the 7T7t scattering 
' lengths [3]. 

ON ! However, the power counting fails when baryons are introduced [4]. The power counting can be restored in heavy 
> baryon chiral perturbation theory (henceforth referred to as HB) [5] where the heavy components of the baryon fields 
are integrated out. An alternative - the infrared regularization scheme (henceforth referred to as IR) - has been 
proposed [6] , based on the ideas in Ref. [7] . This preserves the chiral power counting and has the advantage that it is 
manifestly Lorentz invariant and avoids the voluminous effective Lagrangian of the HB approach. Generally favorable 
results have been obtained with the IR scheme in a number of applications [8]. By suitable approximation of the 
IR expressions the HB formulae can be regained. A recent review of the HB and IR formalisms has been given by 

£ Meifiner [9]. 

The first fit to the pion-nucleon scattering phase shift data using chiral perturbation theory was carried out in Ref. 
[7], henceforth referred to as I, using the HB scheme. In addition to the nucleon and pion fields, the A resonance 
field was included explicitly since the intent was to fit out to energies in the A resonance region. The calculation was 
carried to 0(Q 3 ), where Q is a generic small momentum scale, and it was found possible to obtain a reasonable fit 
up to energies slightly below the A resonance. Subsequently Fettes and Meifiner carried out several HB studies both 
with [10] and without [11,12] explicit inclusion of the A field. The calculation with the A field to 0(Q 3 ) appeared to 
be a little better than the 0(Q 4 ) calculation without it; for further discussion see Ref. [13]. Fettes and Meifiner [14] 
have also studied isospin violation in the 7r — N system, although here we shall focus on the isospin symmetric case. 
All of these calculations were carried out with the HB approach and it is natural to examine the IR method in this 
context. Becher and Leutwyler [15] have used the IR approach in the sub-threshold region, where the HB scheme is 
inappropriate, and found that the IR representation of the scattering amplitude was not sufficiently accurate to allow 
the extrapolation of the experimental data to this region. The purpose of the present work is to study the IR scheme 
in the physical region at 0(Q 3 ) and compare it with the HB approach in order to see whether an improved fit to the 
phase shift data can be obtained. 

The organization of this paper is as follows. In Sec. II we review our notation for the effective Lagrangian and 
discuss the 0(Q 3 ) calculation of the scattering amplitudes in the IR scheme and their reduction to HB form. Formulae 
for the a term and effective vertex couplings are also given. In Sec. Ill our fit to the phase shift data is described, 
first for the HB scheme and then for the IR approach. Our conclusions are presented in Sec. IV. Expressions for the 
IR integrals mentioned in the text are given in the Appendix, together with their HB reductions. 
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II. FORMALISM 



A. Effective Lagrangian 

In the IR or HB schemes the Feynman diagrams for TtN scattering follow Weinberg's power counting rule [16]. 
Define 

u=l + 2L + J2 V i( d i + h n i- 2 ) > (!) 

i 

where L is the number of loops, Vi is the number of vertices of type i characterized by n>i baryon fields and di pion 
derivatives or pion mass factors. Then a given diagram is of leading order Q u , where Q is a small or "soft" momentum 
scale, for example, the pion mass m, the pion momentum or the mass splitting between the A-isobar and the nucleon, 
S = M A - M. 

We briefly recapitulate the effective chiral Lagrangian given in I. Characterising it by di + \rii, the order Q 2 part is 

C 2 = N(Vp + . 9 a7 m 75^ - M)N + i/> (dprfWU) + \m 2 fltv (U + U^- 2) 

+A» A£A£ + h A (a^N + Na^-A^ + h A A^ l5 a^ a , (2) 

where the trace is taken over the isospin matrices. The isotriplet pion fields enter in the SU(2) matrix 

U(x) = C 2 = exp(2m(x)/U) , (3) 

where f n is the pion decay constant and n(x) = tt(x) • |r, with r denoting the Pauli matrices. The axial vector field 
a fl (x) and vector field v t _ l (x) are defined by 

i 11 
a M = -^d^ - £0„ft) = a* = \a^-T = —d^ir - ^ir[ir, 0„7r] + ■ ■ ■ , (4) 

/ 2 



= 2 + =vl= \v,.t = - \l- j [n, <?„*]+■■■, (5) 

both of which contain one derivative. The covariant derivative on the nucleon field V^N = d^N + iv^N. For the A 
resonance we have introduced an isovector field A M = TA p in terms of the standard 2x4 isospin | to | transition 
matrix and the labels a and b in Eq. (2) are isospin indices. The kernel tensor in the A kinetic energy term is 

A» v = -{Hp - M A )g^ + i(^V v + 7^) - + M A )Y , (6) 

suppressing isospin indices, and the covariant derivative is defined by 

Z> M A„ = d^A v + iv^A v — x A„ . (7) 

We have chosen the simplest form for the 7riVA interaction in Eq. (2) since modifications of the type suggested by 
Pascalutsa [17] or modifications of the standard off-shell Z parameter can be absorbed in the other parameters of the 
Lagrangian [18]. 

The order Q 3 and Q 4 parts of the Lagrangian are: 

C 3 = ^NNtr(d^d^U) - ^Nv^N 

+ ^2^ Vu Nti {a»a v ) + jj™ 2 NN tr (U + - 2) + • • • , (8) 



16M 4 
A 6 



+ T7nr7J^7p{^, V v }r a N tr (r a [D"a^, a"]) 



+ 2l ^m 2 {A M -tr[^((7 - rf) T ]N + Ntv[id^(U - U^)r]-A^ + ■■■ . (9) 
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(a) (b) 

FIG. 1. Tree-level diagrams for ttN scattering: (a) contact interactions; (b) nucleon exchange with the cross diagram 
suppressed 



FIG. 2. One-loop A self-energy diagrams. The open box represents the free A propagator. 



Here we have used the definitions 

= - {dp -iVft) ; v liV = d tl v I/ -d v v li + i[v li ,v v ] = -i[a li ,a v ] , (10) 
D^clv = df,a v + ilv^a,,] ; D a v^ = d a v^ + , (11) 

In Eqs. (8) and (9) the ellipsis represents terms that do not contribute to the irN scattering amplitude. 

B. Scattering Amplitudes 

Following the standard notation of Hohler [19] for irN scattering we write the T matrix as 

T ba = (ir b \T\w a ) = T+S ab + \[n,T a ]T- , (12) 

where the isospin symmetric and antisymmetric amplitudes are 

T± = u(p') [A* + \ (4 + /)B ± ] u(p) , (13) 

and u denotes a nucleon spinor. Here, as shown in Fig. 1, q and q' are the cm. momenta of the incoming and 
outgoing pions with isospin labels a and b respectively. The cm. momenta of the incoming and outgoing nucleons are 
labelled p and p' respectively. The amplitudes A ± and are functions of the Mandclstam invariants s = (p + q) 2 , 
t = (q — q') 2 , and u = (p - q') 2 . 

In Fig. 1 we show the tree level Feynman diagrams arising from the contact terms and from one nucleon exchange. 
The vertex in Fig. 1(a) arises from any of the interactions in £ 3 and £ 4 (except for the Ai and A6 terms), as well as 
the Weinberg term —Nj^v^N of Li. The amplitudes for these tree diagrams were given in I. 



1. A Exchange 

When the A appears as an intermediate state for ttN scattering, the tree- level T-matrix diverges at s = M\. In 
I we argued that the power counting should work only for irreducible diagrams. Therefore we summed diagrams 
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containing the one-particle irreducible self-energy insertions shown in Fig. 2 to all orders so as to replace the free 
propagator by the dressed propagator which is finite. To 0(Q 3 ) this gave for the real part of the A propagator in d 
dimensions 



RG„„(fc) = -A R (k 2 )$ + M A )(P 3 /V - 



p 1/2 ) 



(d-2) 
(d - l)M\ 



-Qi + M A )(P^\ V 



(14) 



where the spin projection operators [20,21], denoted by (f^) v , are generalized by replacing factors of 2 and 3 by 
(d — 2) and (d — 1), respectively. In Eq. (14) we have defined 



A fl (fc 2 ) 



A: 2 



n A (fc 2 ) 



[k 2 - Ml - n A (fc 2 )] 2 + M A r A (fc 2 ) ' 



(15) 



(Obviously the bare propagator corresponds to setting II A and T A to zero.) 

The diagrams which contribute to the self energy are shown in Fig. 2, where the open box denotes a free A 
propagator. Only the g^ v part is needed in leading order and the renormalized self energy is 



Y^(k) =S A (fc)-3?S A (fc) 



_d_ 



3?S A (fc) 



(16) 



The integrals which arise are evaluated in the IR scheme, which is briefly discussed in the Appendix. For the 7- 
matrix algebra we evaluate terms whose leading contribution is at 0(Q 3 ) and we discard terms for which the leading 
contribution is of higher order. We include polynomial terms obtained from the product of a l/(d — 4) singularity 
with (d — 4) factors which were were dropped in I. Following this procedure the real part of the A polarization in the 
s- or u-channel is given by 



n A (s) 



h 2 



(s + 2MM A + Ml) l[f (s, M) - iH> (Ml, M) 



(2), 



-2(s - Ml)M A (M A + M)^Xi^(M 2 , M) 



9 -t(2), 



2(s + 9Ml) i[f(s,M A )-i£>(Ml,M A ) 



(2), 



d 



-20(s-Ml)Ml-^i[l\k 2 ,M A ) 



+ 



13(s- 



Ml) 



36s 2 



k 2 =M^ 



(s + 3M A ) 6sm 2 - (s - Ml) 2 



24s 2 ra" 



(17) 



-(2) 

where the integral T\y is specified in the Appendix. By suitably approximating the above expression the HB form 
given in I is obtained (apart from polynomial terms). The width is zero in the w-channel and in the s-channel it is 



irh 2 A 



12M A s 2 (4irU) 2 



[(a - M 2 + to 2 ) 2 - 4sm 2 ] 2 (s + M 2 + 2MM A - m 2 ) 



(18) 



which gives the exact result [19] on shell at s = Ml. 

The A exchange tree-level diagram for nN scattering is pictured in Fig. 3, where the solid box denotes the dressed 
propagator discussed above (a similar notation is used in Fig. 5 below). The corresponding contribution to the 
T-matrix can be obtained from the expression given in I. 

Including also the loop diagrams, discussed below, the real part of the total T-matrix yields the real part of the 
elastic scattering amplitude, f a , by means of the standard partial wave expansion [22]. Here the isospin-spin partial 
wave channels are labelled by a = (I, 21, 2 J) with I the orbital angular momentum, / the total isospin, and J = / ± | 
the total angular momentum. The phase shifts 5 a are then given by 
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(19) 
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FIG. 3. The A exchange diagram with the dressed propagator represented by a solid box. The cross diagram is not shown. 



where \q\ is the magnitude of the c.m. three-momentum. We shall refer to this as the "S'-matrix" approach. 

An alternative, which has been espoused by Fettes and Meifiner [10-12], is the if -matrix approach. This is intro- 
duced by setting 

fa = tth T ^ where K a = tan S a ■ (20) 

l-i\q\K a \q\ 

The calculated amplitude 5i/ a is then assumed to actually be K a so that an infinite f a on resonance corresponds to 
S = i-7T (for further discussion see [23,24]). This allows the free A propagator to be used everywhere. Thus it is not 
necessary to sum the self energy insertions of Fig. 2 and we can expand A# to first order in IIa- The zeroth order 
term gives the tree diagram of Fig. 3, but with the free propagator, and the first order term is treated in the same 
way as the other third order loop diagrams. We will compare this method with the S'-matrix approach. 



2. One-Loop Diagrams 



A set of one-loop diagrams that contribute to the ttN T-matrix at 0(Q 3 ) is shown in Fig. 4. Here we can use 
the free A propagator, denoted by an open box, since no singularities are generated in the T-matrix. We illustrate 
our procedure by discussing the evaluation of Fig. 4(g), which has the same value as Fig. 4(h). Using the standard 
Feynman rules, we find 



Ug 



-«(p') 



i\9a^a - a- 



d J i -^tT c G vtl (p' + *)Tty" 



i 



M 



M + ie' 



C 2 



(21) 



where /j, is the renormalization scale, T l denotes a component of the isospin | to | transition matrix and the subscript 
I denotes that the integral is to be evaluated using infrared regularization. Rather than using the projection operators 
in (14) it is more convenient to write the free A propagator in the form 



G^(fc) - 



# - M A + it 



d- 1 



(d - 1)M A 



(7„fc M - Kin) 



(d-2) 
(d-l)Mi 



(22) 



To 0(Q 3 ) the f terms in the baryon denominators can be dropped. Then the integral involves l v l p in the numerator 
and to this order only the g vp contribution is needed. Carrying out the 7-matrix and isospin algebra, we obtain 



Ug 



-U(P') 



i\ 1 

T . T 



lg 2 A h\ 



7 M 2 {M + M A )(s -M 2 -m 2 + \t- 2Mi) 



l[f (s, t, M A ,M,M)- ^4 (s - Ml - m a ) I | + 32*>L + In ^ 
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u(p) . 



(23) 



(a) 



(g) 



(m) 



\ 1 



/ \ 



/ T 

(b) 



(n) 



* ' - s 

y 



(c) 



(h) (i) 



(o) 



(j) 



> i 



(d) (e) 



JL 




(k) 



) 

/r 

(P) (q) 

\ JL 



(s) 



(t) 



(u) 



(f) 



r 

a) 



(r) 



FIG. 4. A set of one-loop diagrams which contribute at 0(Q 3 ). Crossed diagrams for (d) to (n) are not shown. 
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FIG. 5. Diagrams with one-loop vertices which contribute at 0{Q i ). Crossed diagrams are not shown. The solid circle in 
(a) and (g) refers to A x and \ 6 vertices, respectively. Each diagram implicitly includes its counterpart where the lower vertex 
is dressed. 



The integral is defined in the Appendix, as is the divergent quantity L. In the usual way the terms involving 

(2L + In can be absorbed in the low energy constants so that the result is indcpcndant of the renormalization 

scale. This expression can be reduced to HB form. Here, as in I, we define the average baryon mass, M = \{M + Ma) 

and 77 = (s — M 2 )/(2M). Then evaluating the expressions to leading order and using the reduction of given in 
the Appendix, we obtain 

T 43 - u{ P ')\[t\ r a ] 27 p^ A u)2 [Mi — Mr) + I(2m 2 - t)] [F 2 (rj, -6) §(„ - 5)] u(p), (24) 

recalling that S denotes the delta-nucleon mass difference. This agrees with the result given in I, apart from the 
polynomial term which we include here, but which was subsumed in the low energy constants in I. To these results 
in the s-channel should be added the contribution in the ^-channel with the replacement s — > u and the interchanges 
a <-» b and q <-> —q'. 

Proceeding in similar fashion the real parts of the T-matrix for the diagrams of Fig. 4 may be calculated. Note that 
the diagrams in Fig. 4(m) and (n) involve the nucleon one- loop self-energy which is real for the energies of interest 
here. As with the A, we make on-shell mass and wavefunction counterterm subtractions to obtain the renormalized 
self-energy, namely 

E^ n (fc) = H N {k) - V N (k) -|rSiv(fc) Qt - M) . (25) 



lf=M 



Since these diagrams do not give singular contributions to the T-matrix we do not sum these self-energy insertions. 

We also need to evaluate in similar fashion the diagrams of Fig. 5, where the solid boxes denote the dressed 
propagator discussed in Subsec. II.B.l. The diagrams of Fig. 5(a)-(f) modify the irNN tree vertex in £2, while 
diagrams (g)-(l) similarly modify the 7rA^A tree vertex. 

In a few cases the IR results can be checked against Ref. [15]. The reduction of the IR expressions to HB form 
also provides a check since the results of I should be reproduced (these, in turn, were checked against Mojzis [25] for 
the cases without a A). In the course of checking we found a phase error in one of the HB results in I. This means 
that the expression given in I for diagram 4(i), which is equal to that for diagram 4(j), should be multiplied by a 
factor of 5. This has little impact on the fits presented in I, but the parameters are modified somewhat. Since the IR 
results will not prove to be satisfactory, we shall not list the rather lengthy expressions for all the diagrams; they are 
available on request. 



C. a Term and Effective Couplings 



We may obtain the nucleon a term from the Feynman-Hcllman theorem, 

m <»> 

The nucleon mass receives contributions from the n 2 term in £ 3 and diagrams (m) and (n) of Fig. 4. In the IR scheme 
these yield 



4:K 2 m 2 9g 2 A Mm 2 - ,„,, . 



2h\(M + M A )m 2 [- /li/f2 A/f ^ t (M 2 - M A )(M 2 + 2M|) 



In the HB approximation this reduces to 



(27) 



4« 2 m 2 97rg 2 m 3 4/^m 2 

CT(0) ^ ~ 4(i^F - (i^ [J(,5) - s] ' (28) 

where the integral J is defined in I, see also the Appendix. This agrees with the result of Fettes and Meifiner [10] 
and, modulo a polynomial contribution in m 2 5, with I. 

The irNN vertex up to one-loop order consists of the tree vertex generated from the axial a M term in C 2 and the 
one-loop diagrams shown in the upper part of Fig. 5 - diagrams (a) to (f). As in I, we can calculate the one-loop 
vertex function T a (fc, k', q), where k(k') is the incoming (outgoing) momentum of the nucleon and q — k' — k is the 
momentum transfer. The ttNN coupling for on-shell nucleons is then obtained from 

u{k')T a {k, k', q)u(k) = g^NN{q 2 )u{k') lb T a u{k) . (29) 

At zero momentum transfer we obtain, to 0(Q 2 ), 

, _ 9*nn{0)U , 2 ™ 2A i , 9 A ™ 2 Uh 2 A m 2 M 

(,ttNN = 77 = 1 775- + 777-; — 7^-7 + 



Mg A g A M 2 2(4tt/ t ) 2 27(M a - M)(4tt/ w ) 2 

2 -2m 2 
4m 2 M 2 



n- ,,,2 [(M 2 -M A ) 2 -2m 2 (M 2 + M A )]- „, , 

I n (M 2 ,M) + ^ ^X n (M 2 ,M A ) 



+ ( Tl^ 2) [(M2 ' Mi)2 ' 6w2M2] } 

- W ^^^ ) { [7MMi + 8Ml " M2(2M + 3MA)] ^ (M2 ' MA) 

+ T2^4 I 35MM A + 37M A - 2M 2 (2M + 3M A )][2m 2 M 2 - (M 2 - M|) 2 ] J . (30) 

With the parameters obtained from fits to the irN phase shifts this allows a test of the Goldberger-Treiman relation. 

Similarly, we can calculate the ttNA vertex from the tree-level h A term and the one-loop diagrams shown in the 
lower part of Fig. 5 - diagrams (g) to (1). Here, in the vertex function r fla (k,k' ,q) the label k now refers to the 
incoming A momentum. The 7r7VA coupling is obtained from 

u(k')T^(k, k', q)u»(k) = g,NA(q 2 )u(k')q^T a u^k) , (31) 

where u^{k) is the A spinor. At zero momentum transfer we obtain, to 0(Q 2 ), 

9*na(0)U_, 2m 2 A 6 4ff 2 i M 2 - (3) 2 ^ ^ 65g A h A (M + M A )r 



X$>(Mi,M,M) 



Mh A h A M 2 (4ir,U) 2 ' ; m=o 108M(4 7 r/ 7r ) 2 

2h 2 A M(M + M A ) f (3) 2 1 [(M 2 -M 2 ) 2 -6m 2 M 2 



9(4ttA) 2 I 
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50h A M A (M + M A ) (3) 2 

81(47r/ff)2 < X 12 (M A , M A , Ma) 



m=0 



(15M + M A )[(Mj - M 2 ) 2 - 6m 2 M 2 ] ] 

360M 4 M A J 1 J 

Note that for the first two integrals 2^ 2 both real and imaginary parts need to be considered. The value of g W NA 
is complex because the intermediate pion and nucleon states for Fig. 5(i) and (j) can go on shell. We therefore 
define CttNA = |5irJVA(0)/ ff /(M/iA)|- The HB reduction of Eqs. (30) and (32) agrees with the expressions in I up to 
polynomial terms. 



III. RESULTS 



As fixed input parameters, we use the standard baryon and pion masses: M — 939 MeV, M A = 1232 MeV, and 
m = 139 MeV. We also take [26] /„• — 92.4 MeV from charged pion decay, <7a = 1-26 from neutron (3 decay, and 
h A = 1-46 from Eq. (18) for T A (M A ) using the central value of the A width, T A = 120 ± 5 MeV. We have ten low 
energy constants: (3 n , k v , k\, k 2 , and Ai to A 6 , plus the 7rAA coupling, Jia- These are obtained by optimizing the fit 
of our calculated irN S- and P-wave phase shifts to the 2001 data of the VPI/GW group [27]. Since errors are not 
given we compared results obtained by assigning all the data points the same relative weight to those obtained using 
weightings suggested in the literature [10,28]. The differences were small so we used the same relative weightings, 
minimizing 

* N / ^theory ^-experiment \ 2 

y 2 = ir Si Zil (33) 

A J\[ / / I q Qg^cxperiment I ' V > 

where = sin 25, for the S- matrix calculation and £j = 2 tan Si for the K- matrix calculation. 



A. HB Results 



Our fits to the S- and P-wave phase shifts in the HB calculation are shown in Figs. 6 and 7 using the S'-matrix 
and if-matrix methods, respectively. Here the solid curves correspond to an unconstrained fit. As in I and Ref. [10] 
it was possible to fit the phase shift data out to a pion cm. kinetic energy, e, of 100 MeV, slightly below the A 
resonance at e = 127 MeV. More precisely the data were fitted from ^fs = 1090 to 1200 MeV in 3 MeV steps. The 
short horizontal section in Fig. 6 indicates a small region where unitarity, which is only enforced perturbatively, is 
slightly violated and a phase shift cannot be determined; this, of course, does not occur in the iVmatrix approach, 
by definition. The corresponding parameters are collected in Table I where, in common with previous work, it is 
observed that the deduced values of the a term are much larger than the oft-quoted figure of 45 ± 8 due to Gasser 
ct al. [29]. Therefore we performed two additional fits where k 2 was constrained to produce values of the a term of 
75 MeV and 45 MeV. These are indicated, respectively, by the dotted and dashed curves in Figs. 6 and 7. For the 
if -matrix case the dotted curve is barely distinuishable from the solid curve since the change in the a term is small 
here. 

As expected the S-matrix fits in Fig. 6 are quite similar to those given in I. The fitted parameters, however, differ 
since here we include the polynomial terms and we subsume the rcnormalization scale dependance into the low energy 
constants. In many of the partial waves the K- matrix fits (Fig. 7) are a little better, particularly for <5p33, however 
Spsi begins to deviate significantly from the data at e ~ 80 MeV. Something of the same trend is visible in the work 
of Fettes and Meifiner [10]. The 0(Q 2 ) low energy constants given there are related to our constants by c\ = K2/M, 
C2 = ki/(2M), c 3 = 2(3 7T /M and c 4 = K n /M. Their fit l 1 * to Matsinos' data [30] produces values which are fairly 
similar to ours when the different value of the ttNA coupling, Ha, and the constraint <r(0) = 59 MeV are taken into 
account. (We should note, however, that substantial differences are seen for their fit to the Karlsruhe data [31].) The 
0(Q 2 ) constants in Table I are reasonably similar in the S- and iVmatrix calculations, however sizeable differences 
appear at C(Q 3 ). Further they are not of natural size (unity). These features may be due to the strong cancellations 
which occur among the various terms and the fact that the parameters are not uniquely determined since it is possible 
to obtain similar values of x 2 when <r(0) is constrained. A significantly larger x 2 is only obtained for the S'-matrix 
case with <r(0) = 45 MeV. 
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FIG. 6. The S- and P-wave phase shifts using the S-matrix method as a function of the pion cm. kinetic energy, e. The 
solid curve corresponds to an unconstrained nucleon a term, while for the dotted and dashed curves it is constrained to be 75 
MeV and 45 MeV, respectively. The data are from Ref. [27] 



10 




e (MeV) e (MeV) 



FIG. 7. As for Fig. 6, but using the 7f-matrix method. 
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TABLE I. Low energy constants and effective couplings obtained from the S- and /^-matrix fits with unconstrained and 
constrained values of the nucleon a term. 





S'-matrix 


if-matrix j 


ft 


4.78 


4.46 


4.13 


4.06 


4.06 


4.05 


Kir 


-10.1 


-9.02 


-7.76 


-7.96 


-7.79 


-7.21 


Kl 


-10.6 


-9.80 


-8.82 


-9.36 


-9.29 


-9.04 


K2 


—0.453 


—0.124* 


0.241* 


—0.234 


—0.124* 


0.240* 


Ai 


23.7 


16.6 


8.20 


9.11 


7.77 


3.13 


A 2 


-20.1 


-18.0 


-15.9 


-14.9 


-14.7 


-14.4 


A 3 


-19.5 


-16.8 


-13.4 


-14.1 


-13.5 


-11.6 


A 4 


-21.0 


-19.0 


-17.1 


-16.9 


-16.8 


-16.5 


A 6 


9.14 


8.03 


7.20 


6.03 


6.00 


6.01 


A 6 


0.661 


-3.01 


-7.08 


-7.32 


-8.03 


-10.4 


flA 


1.24 


0.918 


0.543 


0.612 


0.548 


0.329 


x 2 


0.26 


0.38 


0.77 


0.92 


0.94 


1.2 


<t(0) (MeV) 


102.1 


75.0 


45.0 


84.0 


75.0 


45.0 


CttNN 


1.088 


1.056 


1.021 


1.050 


1.041 


1.012 


CitJVA 


0.905 


1.024 


1.156 


1.160 


1.183 


1.260 



* Constrained. 



TABLE 11. S'-wave scattering lengths and P-wave scattering volumes in units of m 1 and m 3 , respectively 





S-matrix 


if-matrix 


Experiment [ 


(7(0) (MeV) 


102.1 


75.0 


45.0 


84.0 


75.0 


45.0 




60 


0.0051 


-0.0061 


-0.0180 


0.0016 


-0.0024 


-0.0155 


-0.0012 ±0.0010 [35] 
-0.0034 ± 0.0007 [36] 


61 


-0.0830 


-0.0826 


-0.0824 


-0.0847 


-0.0846 


-0.0843 


-0.0895 ±0.0016 [35] 
-0.0918 ±0.0013 [36] 


an 


-0.0753 


-0.0754 


-0.0752 


-0.0769 


-0.0766 


-0.0757 


-0.078 ± 0.002 [31] 


ai3 


-0.0290 


-0.0289 


-0.0282 


-0.0281 


-0.0280 


-0.0276 


-0.030 ± 0.002 [31] 


d31 


-0.0392 


-0.0398 


-0.0400 


-0.0378 


-0.0379 


-0.0382 


-0.046 ± 0.007 [37] 


a33 


0.2180 


0.1984 


0.1787 


0.1966 


0.1919 


0.1766 


0.205 ± 0.004 [37] 



The ratios of the effective irNN coupling constant to the bare value, C^nn, in Table I show a Goldberger-Treiman dis- 
crepancy of a few percent. Schroder et al. [32] have recently determined a precise result, namely CttAW = 1-027+ 008 . 
A similar value is obtained from the latest ttNN coupling constant obtained by the George Washington Univer- 
sity/TRIUMF group [33]. These values favor a a term of between 75 and 45 MeV which is in accord with several 
recent analyses. In Ref. [33] 64 MeV is quoted with an error of about 10%. A sum rule determination by Olsson 
[34] gave 55 MeV with a 16% error. Finally Schroder et al. [32] have indicated that the value of 45 MeV extracted 
[29] from the Karlsruhe data should be increased by 13 MeV giving 58 MeV. Lastly we note that since remains 
fairly close to unity the coupling is not changed too much from the bare value, which seems intuitively reasonable. 

It is also interesting to examine the threshold results. We give in Table II the S'-wave isoscalar and isovector 
scattering lengths, 60 = ( a i ± 2a 3 )/3 and 61 = (a 3 — ai)/3 (with the notation a 2 i) and the P-wave scattering volumes 
(fl2i 2j)- The experimental values come from various sources. The S'-wave isoscalar scattering lengths were obtained 
[35,36] from the pionic atom data [32] with an improved treatment of the 7r-<i scattering length. The I = \ P-wave 
scattering volumes are the old Karlsruhe values [31], while the I = § volumes are from a recent analysis by Fcttcs and 
Matsinos [37]. Our calculated results for bo at the extreme values of <r(0) are clearly not favored by the data which 
prefer a value in the middle. Notice that, at this level of accuracy, there is a noticeable difference between the S- 
and if -matrix results. The value of 1233 appears to be somewhat low for the cases where cr(0) = 45 MeV. Apart from 
this, the remaining P-wave results and the values of b\ show little sensitivity to the calculation employed and are in 
reasonable agreement with the data. Of course this is not surprising since the low energy phase shifts are included in 
the fit. 
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TABLE III. Low energy constants, a term and effective couplings obtained by fitting to e = 40 MeV using the S- and 
if-matrix in the IR scheme and the if-matrix in the HB scheme 







IR 




HB 




S matrix 




t\ — m a t~ T"i v 

1 i 111CLL1 IjV 


l\ — m a tri y 

1 V 111CLL1 IjV 


a 


4.23 




O.O-L 


6.12 




— 10.4 




— 12.9 


— 12.0 


K-\ 

rvi 


— 10.4 




— 15.1 


— 14.5 


K2 


-0.121 




-0.244 


-0.287 


Ai 


22.0 




33.6 


31.5 


A 2 


-24.4 




-26.9 


-23.7 


A3 


-38.2 




-45.3 


-23.9 


A 4 


-16.9 




-18.9 


-24.4 


A 5 


18.0 




18.8 


7.74 


A 6 


19.6 




19.1 


-6.00 


h A 


0.500 




0.664 


1.50 


x 2 


0.075 




0.066 


0.025 


<r(0)(MeV) 


108.9 




119.0 


88.4 


CirJViV 


0.790 




0.699 


1.043 


CirJVA 


0.378 




0.381 


1.079 



TABLE IV. Comparison of IR results for the real parts of various integrals and IIa in the s-channel with the corresponding 
HB approximation. 





y/i = 


1200 MeV 


>/i = 


1100 MeV | 


Integral 


IR 


HB 


IR 


HB 


In (s, M A ) 


0.338 


0.357 


0.148 


0.120 


i^>{s,M A )/M 2 


0.00193 


0.00199 


0.000314 


0.000192 


T^(s,M A ) 


-0.00127 


-0.00413 


-0.00797 


-0.0103 


M 2 T^(s,Q,M a ,M,Ma) 


0.0798 


0.0453 


0.1020 


0.0568 


n A (s)/M 2 for h\ =0.1 


0.000287 


0.00546 


0.0230 


0.0696 



B. IR Results 

Turning to the IR scheme the results are disappointing. It was only found possible to fit the phase shift data in 
the low energy regime out to e = 40 MeV. Specifically we fitted for y/s from 1090 to 1123 MeV in steps of 3 MeV. 
The IR 5*-matrix results are shown by the dashed curves in Fig. 8; here the restriction Jia < 0.5 was imposed, see 
the discussion below. The if-matrix results are denoted by dotted curves in the figure. For contrast the solid line 
shows the HB results for a fit over the same energy range; the HB results from the S- and if -matrix methods are 
very similar and we choose to display the If -matrix results. While the HB phase shifts follow the data approximately 
for e > 40 MeV, in a number of cases the IR results rapidly diverge from the experimental values. The parameters 
corresponding to these results are collected in Table III. While the HB parameters differ from those given in Table 
I, particularly those of 0(Q 3 ), we caution that they are not well determined by a fit over this limited energy range. 
Notice that, for the if -matrix cases, x 2 is substantially larger in the IR scheme than the HB one. In the IR case (-knn 
differs substantially from unity, implying a huge Goldberger-Treiman discrepancy. This is unphysical and is a signal 
of the contortions that the parameters are going through in order to fit over even this limited range of energies. The 
quantity CttJVA is also much less than unity implying an effective nNA coupling very different from the bare value. 

We have, of course, carried out a number of additional IR calculations with a view to improving the results. For 
example constraining the value of the a term to be 75 or 45 MeV increases the value of x 2 > but the results are 
qualitatively very similar to those discussed above. In order to see whether the problem was due to the A resonance 
contributions we removed them entirely by setting the coupling Ha to zero. As expected the parameters changed 
substantially, but the fit to the phase shifts remained qualitatively the same and in fact x 2 increased somewhat. Thus 
we have not been able to qualitatively improve the IR results discussed in detail above. 

In order to further investigate the reason for the problem with the IR approach we compare the IR polarization Ha 
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FIG. 8. The S- and P-wave phase shifts from the S-matrix in the IR scheme (dashed curve) and from the if-matrix in the 
IR scheme (dotted) and HB scheme (solid) as a function of the pion cm. kinetic energy, e. The data are from Ref. [27] 
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TABLE V. Comparison of IR results for the real parts of various integrals and ITa in the tt-channel, for u = u min , with the 
corresponding HB approximation. 



Integral 


= 1200 MeV 

ylimin — OUl./ IVleV 

IR HB 


yfs = 1100 MeV 
V^min — <oy.o iviev 
IR HB 


Tn (w, Ma) 
l[f(u,M A )/M 2 

M 2 T[f(u,0,M A ,M,M A ) 
n A (u)/M 2 for h\ = 0.1 


-6.51 -1.42 
1.34 0.0826 
0.840 0.110 
0.597 0.167 
-0.751 0.319 


-2.39 -0.934 
0.173 0.0316 
0.180 0.0500 
0.321 0.130 

0.0213 0.262 



of Eq. (17) and a few IR integrals (denned in the Appendix) with their HB counterparts. In Table IV we make the 
comparison in the s-channel. For the most part the IR and HB results are quite comparable. The picture changes for 
the tt-channel in Table V, where for a given s we have chosen the minimum value, tt m ; n = 2(M 2 + to 2 ) — s, so as to 
enhance the contrast. Here the IR integrals are much larger than the corresponding HB results, by a factor of 6 on 
average. Further Ua(u) changes sign. Since A^(tt) = l/[tt — M A — ILa('u)] and tt — M A is negative this change of sign 
can produce a pole; this was the reason for our restriction of Ha < 0.5 in the S'-matrix IR fit. Of course the problem 
does not arise if n A (tt) is expanded out of the denominator as in the if -matrix approach, although the mathematical 
justification for doing so is slight. The problem does not arise in the HB case where IIa(u) is always positive. 

It is also useful to make the point graphically by studying the basic integral involving one baryon and one meson 
propagator. Figure 9 shows Xn(s, Ma) as a function of y/s. It is observed that the IR integral (solid line) and the HB 
approximation thereto (dashed curve) agree quite closely. Figure 10 shows the corresponding result in the tt-channel; 
Zii(«min, Ma) as a function of y^s. It is seen that the magnitude of the IR integral (solid line) is much larger than 
that of its HB counterpart (dashed curve); note the ordinate scale here in comparison to Fig. 9. The IR integral is 
governed by the value of £ (see Appendix) which in the s-channel is (s — M^ + m )/(2m^/s). In the HB approximation 
this is replaced by (s - M A )/(2mM) = (r? - 5)/m, in the notation of I. At y/s = 1200 MeV, x is -0.18 (-0.26) in 
the IR (HB) case. Correspondingly in the tt-channel, at ^/zt m in = 601.7 MeV, x is —6.79 (—3.83) in the IR (HB) case 
and this large difference is the heart of the problem with the IR approach. However, it is disquieting that the rather 
reasonable HB results are achieved by the replacement of y/u in the denominator by the square root of average baryon 
mass, VM, since this can be quite a rough approximation. 

IV. CONCLUSIONS 

We have performed a least-squares fit to the S- and P-wave phase shift data using chiral perturbation theory to 
0(Q 3 ) with explicit tt, N and A fields. The heavy baryon formulation used here represents an improvement compared 
to I in that the polynomial terms were included explicitly and the renormalization scale dependance was subsumed 
into the low energy constants. We contrasted an "S'-matrix" method (as in I) with a if-matrix analysis (as in Ref. 
[10]). The former seems to be more physical since it employs dressed A propagators which are finite, while the latter 
uses bare propagators which give a divergent scattering amplitude at the A resonance energy. Broadly speaking 
the two methods give similar results and allow a fit to the phase shift up to a cm. pion kinetic energy, e, of 100 
MeV. In detail, however, there are some differences. For instance, the value of the a term (if unconstrained) and the 
Goldberger-Treiman discrepancy are smaller with the if -matrix. 

Our main interest was to see whether the infrared regularization scheme, which is manifestly Lorentz invariant, 
could provide an equally satisfactory account of the phase shift data. It could not. We were only able to fit the data 
up to e = 40 MeV and the resulting Goldberger-Treiman discrepancy was ridiculously large. The difficulty was not 
due to the explicit inclusion of the A resonance, but rather stemmed from the tt channel where the magnitudes of the 
integrals were much larger than the corresponding heavy baryon ones. This leads to uneasiness regarding the latter 
since they are obtained by approximating the former. The heavy baryon scheme thus pushes significant contributions 
into higher orders where they are combined with the other terms of that order. The only information we have as 
to whether the net effect is small comes from the comparison of third and fourth orders carried out by Fettes and 
Meifiner [12]. In some cases there are sizeable differences, although these authors suggest that in most cases the chiral 
series appears convergent. Overall, in our opinion, the description of pion-nucleon scattering via chiral perturbation 
theory is not yet in a satisfactory state. 
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APPENDIX 



Integrals involving only meson propagators need no discussion since the standard results are used and integrals 
involving only baryon propagators give no contribution in the IR or HB schemes. Infrared regularization is needed 
for integrals involving both baryon and meson propagators. Here we give only those integrals that are referred to 
in the main text, together with their HB reduction. The integrals needed can be derived from the basic infrared 
integral which has been discussed by Becher and Leutwyler [6] (note that our conventions differ from theirs), see also 
I. Explicitly this is 



tjj, 



A-d 



d d £ 



1 



(2w) d [(£ + k) 2 -M 2 + ie][£ 2 -m 2 + ie] 
Ida i 



= ijj, 



d f dH 
J (2ir) d 



(2ir) d [2k ■ £ + k 2 - M 2 + m 2 + ie] [f 2 - m 2 + ie] 
I^ 1 (s,M) + (4tt)- 2 T 11 ( S ,M) , 



(Al) 



where s = k 2 . The total integral, without the subscript /, will have, in addition to the infrared part, a regular part 
but this can be discarded since it can be expanded with the various terms absorbed into the low energy constants. 
The divergent IR contribution is 



where L 



(s - M 2 + m 2 ) 



2s 



2L 



1 



1 



+ 7-1 



32tt 2 

with 7 denoting Euler's constant. The real part of the finite piece is 

-x + 2^1 



(4tt) 2 
In 4-7T 



2 

m 

In -^r 



(A2) 



3?X n (s,M) 



m 

■5' 



\-x)] 



+ Vx 2 - 1 In 



M' 



x-{-\/x 2 — 1 



x — y/x 2 — 1 



\x\ < 1) 

M > i) 



2M 



(A3) 



where x = (s — A4 2 + m 2 ) / \2m^/s) . The last line here makes contact with the HB integral J in I for which s/s in the 
denominator is approximated in leading order by the average baryon mass M = \{M + Ma). The imaginary part of 
the IR integral is 



--V(s -M 2 + m 2 ) 2 - Asm 2 6{s - [M + mf) 



(A4) 



where denotes the Heaviside step function. 

In the following we shall implictly give only the real part of the IR integral and the ie in the integrand denominators 
will be suppressed. For the integral 



IjJL 



J I (27T 



) d [(£ + k) 2 - M 2 ] [£ 2 - m 2 ] 



= l{ 1 1 \s,M)k fl 



The finite part is (4 7 r)- 2 xj 1 ) , with 



T[\\s,M) = - 



s - M 2 + m 2 



2s 



Tu(s,M) 



M 



2 1 



2M 3 



2M 



(A5) 



(A6) 



with the last expression giving the HB reduction. 
For the integral 
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IjJL 



d d P 



P P 



It {2w) d [(£ + k) 2 - M 2 ][P 2 - to 2 ] 
The finite part of the first term is (A-k)~ 2 T^ , with 

"4m 2 s - (s - M 2 + to 2 ) 2 



(3), 



(A7) 



T[f(s,M) 



12s 



M) {S ^f ] [6m 2 s -(s- M 2 ) 2 ] 



1 



2 (s-M 2 ) 2 



AM 2 



J 



M' 



2M 



36s 2 
(s~M 2 ) 
18M 2 



*»»-<- ^ 



2M 2 



(A8) 



with the last line giving the HB reduction. 

One of the integrals involving one pion and two baryon denominators is 



. 4 -d f d d P 

* M J j (2ir) d 1(1 + p< + q>) 2 - 



f P 



M 2 A ][(P + p') 2 - M 2 B ][P 2 - m 2 



= l[l\s,MA,M B )g^ + ■■■ 
The finite part is (4ir)- 2 i[f , with 

( d 



(A9) 



T[f(s, M A , Mb) 



\dM 



dM% 



dzl[f(k 2 ,M) 



where k 2 = m 2 z 2 + (s - M 2 - m 2 )z + M 2 , 
and M 2 = m 2 z 2 + {M A -M 2 B -m 2 )z + M B . 



(A10) 



This was evaluated by performing the differentiation analytically and carrying out the integration numerically using 
Gaussian quadrature. It is reduced to HB form by dropping terms involving m 2 from the expression k 2 — M 2 + m 2 
and setting k 2 ~ M 2 in denominators. This allows the integral to be performed, giving 



T[f(s,M A ,M) 



1 



6M 2 
1 



[F 3 ( V -6,0)-m 2 + 2 -(r 1 -S) 2 ] , 



where F 3 {a, 0) = — - [(a 2 - m 2 )J(a) - (0 1 - to 2 ) J(/3)] 

and we recall the definitions of I: 2Mr\ = s — M 2 and 5 = Ma — M. 
The integral involving one pion and three baryon denominators is 

P P 



(All) 



. 4 _ d f _d^_ 

J I (2nr[(e + p' 



) 2 - M\] [(P + P ' + q') 2 - M%] [(P + p) 2 - M 2 C ] [P 2 - m 2 ] 



= l{f(s,t,MA,M B ,Mc)g^ + ■■■ ■ 
The finite part is (An)- 2 i[f , with 

d d 



(A12) 



iil\s,t,MA,M B ,Mc 



dM 2 A 



dM% 



dM 2 n 



dxx I dyl {2 ^(k 2 ,M) 



where k 2 = [s - M 2 - to 2 + to 2 :e(1 - y)] x(l - y) - txy(l - x) + M 2 
and M 2 = m 2 x 2 {\ - yf + [M A - M% + to 2 - t(l - a;)] xy 
+ (M 2 B - M 2 C - m 2 )x + M 2 C . 



(A13) 



This was evaluated with a combination of numerical differentiation and Gaussian numerical integration. The HB 
reduction is carried out by making similar approximations to those mentioned above and this allows both integrals to 
be performed. For the cases mentioned in the text this gives 
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and 



T[f(s, t, Ma, M, M A ) -_ [F, („, -5) + §(„ - 25)] 



r( 3 )c e + M. M M, \ L 

12M 3 

(a 2 - m 2 ) J(a) + (2/? 2 - 3a/3 + m 2 ) J(/3) + (a - /3)(2/? 2 - m 2 ) 
where Fi («,/?) = — - ^ 2 , (A14) 



1 

12M 3 

(a 2 - m 2 )/3J(a) - ((3 2 - m 2 )aJ(/3) + 7rm 3 (a - j3) 



Tg\s, t, Ma, M, M) -> [F 2 (r?, -<5) + §(??- *)] 



where F 2 (a,/3) = ^ " K ' ^ i ^. (A15) 

a(3{a - (3) 
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